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Introduction 


This is the sixth quarterly report on Contract NAS5-27382 entitled, 
"Spectroradlometric Calibration of the Theaatic Mapper and the 
Multlspectral Scanner System.” In this report, we summarize the 
experiments performed during the April 1984 trip to White Sands, New 
Mexico. This includes the field testing of the newly built 
spectropolarlmeter, both in the solar radiometer mode and helicopter 
flight test. Also included is a preprint of a paper entitled "Effective 
Bandwidths for Landsat-4 and Landsat-D' Multlspectral Scanner and 
Thematic Mapper Subsystems". This is to be published in IEEE Transactions 
on Geoscience and Remote Sensing , May 1984. 


Field Measurements 

As Thursday, April 19, 1984 was the first overpass of Landsat 5 in 
which Imagery was to be recorded over White Sands, we went to the area 
to take calibration measurements. Those participating were Barbara 
Capron, Ken Castle, Ron Holm, Ray Jackson, Carol Kastner, Jim Palmer, Amy 
Phillips, Richard Savage, and Phil Slater. It was cloudy and quite windy 
on the day of the Landsat 5 overpass and we were not able to accomplish 
all our objectives. Neither a complete set of solar radiometry 
measurements, nor reflectance measurements of the gypsum were made on 
April 19. Thus no data are as yet available for ground based calibration 
of this Thematic Mapper (TM) sensor, and In addition Langley plot data 
from the radiometer and Castle spectropolarlmeter could not be compared. 
Three major goals, however, were accomplished, and are summarized below. 
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Castle Spectropolarlaetera 

This was the first opportunity to test the newly built Castle 
spectropolarlmeters In the field. On Thursday one of the Instruments 
was used in the solar radlometry mode. As the sun came in and out from 
behind the clouds the neutral density filters automatically switched 
Into place. The instrument cycled through each of the 10 narrow band 
filters, recording the time of day at the start and finish of each cycle. 
There seemed to be a minor problem decoding one of the neutral density 
filter positions, which is now being corrected. 

The second of the two instruments was flown, also on Thursday, 
onboard an Army helicopter. By measuring the radiance incident on this 
sensor, it is hoped that the Herman code output can be verified for an 
intermediate altitude of 10,000 feet above sea level. In order to mount 
the radiometer on the helicopter, a U-shaped bracket was built at the 
heliport hanger. The bracket had slits in the two side walls which the 
radiometer's mounting handles slid into. A slot at the bottom allowed 
the nose of the radiometer to view straight down. In addition the 
instrument could be tilted to any angle to take measurements, and also 
rotated 180' in place to check Instrument settings while in flight. 

On Thursday the helicopter departed the Holloman AFB helipad at 
9:33 a.m. from a ground altitude of 4110 ft. Upon arriving at Chuck site 
we landed momentarily to remove the helicopter windows. Departing 
again at 9:55 a.m. we climbed to 10,000 ft altitude by 10:01. Just as we 
were in position to take readings, it was discovered that the radiometer 
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was frozen In a 45' position (the crew chief had tightened it down too 
hard just before takeoff). There was nothing that could be done but fly 
east of the test sight, so that the instrument appeared to be aimed in 
the right direction (no precise alignment technique was available). Ken 
Castle then cycled through the ten spectral filters on his instrument. 
The helicopter was repositioned, and a second data set was taken. During 
these two measurements the helicopter flew between 9,900 ft and 10,200 
ft (a thermal prevented us from stabilizing to a more precise altitude). 

A third data set was taken while the helicopter was descending to the 
6,000 ft level. A fourth data set was initiated, only to find that the 
instrument was not shifting filters. The entire flight lasted 1.5 hours. 
Upon landing the instrument was inspected. A loose screw was removed, 
after which the instrument appeared to be totally functional. 

Overall the instrument performance was satisfactory. Washers were 
added to the mounting bracket to prevent the instrument from freezing in 
place during the next flight. The following describes their performance 
in the solar radiometer and helicopter modes in more detail. 
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Solar Radiometer 

The instrument was operated with a five degree field of view in 
order to make it easier to track the sun. A tripod was used as a 
support but it was found to be very difficult to position the 20 lb 
instrument and maintain its direction of viewing. (A much sturdier and 
sophisticated alt-azimuth mount is currently being constructed.) The 
instrument also became locked, due to a software error, into a 
configuration where it was always looking through an ND-4 filter (the 
highest density filter on board). The results below show that the sun 
does not saturate the system in this configuration implying that the 
pre-design calculations were valid. A sample data run is shown below. 
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where: pol represents the polarization state transmitted for each 

measurement 

OCN is the digital count output of the spectropolarlmeter 
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These data are a factor of two less than the saturation level for 
the instrument, thus providing a reasonable margin of comfort. 


Helicopter performance 

The primary test proved the Instrument would operate successfully 
when hard-mounted to a helicopter in flight. The second test was to see 
whether it would be able to detect the ground radiance at a height of 
2000 m above ground level (AGL). The instrument was configured to have 
a 1 degree field of view in order to give a ground resolution of around 
30 m. The following table shows the predicted values and a typical set 
of the measured values (the two orthogonal polarization measurements 
were added together to get a total value). Due to the problem with the 
heliocopter mount mentioned earlier, the instrument was only able to 
point at a 45 degree angle to the ground. (Note: the ND«0 filter was 
used for this test case.) 
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where: E is the exo-atmospheric lrradiance 

T are typical optical depths 
p are typical reflectances of gypsum 
L is the ground radiance (E exp(-TsecQ)/ it 


p det 18 the power on Che detector and is calculated 
by multiplying L by the number 3.086E-04 which is a 
correction factor incorporating the effects of 
transmission, solid angle, and aperture 

R is the detector responsivlty 

I s ig is the current leaving the detector 

DCN P re( j is the expected digital counts out of the 
instrument 

DCN A ct is the actual digital counts measured 

As can be seen, there is very good agreement between experiment and 
theory for this simple test case. Also, all the actual values are 
greater than 100 so that the S/N ratio is better than 100:1 for less 
than a 11 noise error. In general, the instrument will be operating 
with a 15 degree field of view when making ground observations so it 
should be able to see surfaces with lower reflectances. The most 
notable observation concerning this aspect of the instrument operation 
was that it was able to successfully hold the respective filter, 
aperture and polarizer positions while being vibrated by the helicopter 
platform. If it can do that, then it should operate very well in a more 
benign environment. 

Additional data have been collected with the instrument configured 
with a 5 degree aperture and a full set of polarization measurements 
made over all the wavelengths. The instrument was pointed at right 
angles to the sun near sunset and the data show the sky to be 40% 
polarized in the blue and dropping to 20 X polarized in the near IR. All 
the measurements appear to be consistent and so the instrument is 


- 7 - 


actually functional aa a spectropolarlmeter. These laat neaaurementa 
were taken here in Tucson. 

Air photography 

The first helicopter flight of the week was taken on Tuesday. It 
had the dual purpose of photographing the test site to be used on 
Thursday, and inspecting the entire area by air to determine other 
potential sites. The flight departed JFK heliport at 10:45 a.m., but was 
forced to return 10 minutes later when it was discovered that there was 
a scheduling conflict for the air space we had requested. After being 
rescheduled a second attempt was made. This time we were airborne 
from 11:10 a.m. to 12:06 p.m. The first area we investigated was Parker 
site. This is a circular area, at least 200 m in diameter, which is 
graded flat and used as an impact area. It was brownish in color, with 
different hues at different radii. Because of the nonuniformity, it is 
unsuitable as a test area. The pilot mentioned that other slmiliar, but 
larger, impact areas existed further north (close to the lava beds). We 
plan to investigate them on a later flight. 

A quick trip was then made to the Northrup strip /Cherry site area. 
In general the sand area looked quite mottled, with well defined regions 
of gray intermixed with highly reflecting regions. Cherry site looked 
quite uniform. In comparison to many other regions. The photographer 
took 2 rolls of positive (slide), and two rolls of negative film, ASA 160 
and 180 respectively. All pictures were shot at a focal length of 70 
mm, 1/500 sec, and between F/8-11. The pictures were processed that 
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same day, after which It was discovered that they were overexposed and 
contained little useful Information. 

On Thursday the same films, focal length and shutter speed were 
used. This time they were shot at F/22-16 for the slides, and F/ 16—1 1 
for the negative film. Photos were taken of Chuck site at both 10,000 ft 
and 6,000 ft elevation. Our sites to either side of the road looked 
quite uniform. The corners of our two 4X4 pixel area were well marked, 
allowing the photographer to zoom in on this area. 

Test Site Selection and Lay Out 

Our only other opportunity to do a calibration was on January 3, 
1983, using an Image from the TM onboard Landsat 4. There was some 
difficulty in identifing the location of Cherry site. Because of this it 
was decided to define a new test area, one which could be accurately 
located on the Landsat image. We therefore changed our test area from 
Cherry site to Chuck site. This latter area was 3.4 miles to the east of 
Cherry site. The new site was selected because of Its uniformity and 
because it was located at the intersection of a nearly 100' bend in the 
road. As it was located along the same straight road as Cherry site, the 
two 4X4 pixel areas were aligned at an angle to the road from 
computations made for Cherry site, after inspecting the previous image. 
The road is known to lie 32*N of east, hence the pixel edge was adjusted 
such that it fell on a true east/west line. The new 4X4 pixel areas 
should now be closely aligned with the TM scan direction. 


We would Ilk* Co thank tha photographer, Frank Trevino, and the 
pilot* Jack Raea and Lt. Kayaa for their co-oparatlon and assistance. 
Special thanka la due Richard Savage and the management of ASL for the 
time they epent in helping ua organize our trip and making aura our vlalt 
waa a* productive as possible. 
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EFFECTIVE BrfWDWIDTHS FOR LAND SAT -4 AND LANDSAT-D' MULTI SPECTRAL 
SCANNER AND THEMATIC httPPER SUBSYSTEMS 

Jamas M. Palmar 
Uni varsity of Arizona 
Optical Sciancas Cantar 
Tucson, AZ 85721 

ABSTRACT 

Tha spectral bands of the Mu 1 1 i spec tral Scanner and Thematic 
Mapper subsystems of Landsat-4 and Landsat-D' have bean analyzed 
using a bandwidth normalization technique based on analysis of 
the moments of the spectral responsitivy curves. The results 
include the effective wavelength, the bandpass, the wavelength 
limits and the normalized responsivi ty for each spectral channel. 
In addition, temperature coefficients for TM PF Channel 6 have 
been derived. The moments normalization method employed yields 
sensor parameters whose derivation is independent of source 
character isterics (i.e., incident eolar spectral irradiance, 
atmospheric transmittance or ground ref 1 ectance) . The errors 
expected using these parameters are lower that those expected 
using other normalization methods. 


INTRODUCTION 


Relative spectral responsivi ties for the Mul tispectrat Scanner 
(MSS) and th# Thematic Mapper <TM> subsystems of Landsat-4 and 
Landsat-D' were obtained from Markham and Barker ( 1983a) and 
Markham and Barker (1983b), respectively. These data, provided in 
tabular form vs. wavelength, were analyzed using a recently 
described normalization technique (Palmer and Tomasko, 1980). 
This technique, based upon an analysis of the moments of the 
spectral responsivity curves, yields effective center wave- 
lengths, bandpasses, passband wavelength limits and equivalent 
squareband responsivi ties. The'*? parameters, when applied to the 
radiometric analysis of continuous sources, give results that ar-e 
more accurate than analyses using conventional normalization 
methods. Each of the 24 detectors of the protoflight (PF, flown 
on Landsat-4> and flight (F, flown on Landsat-D') MSS channels 
has been analyzed, and means and standard deviations have been 
computed for each band (6 detectors) . Comparisons are shown 
between the previously published bandwidths and center wave- 
lengths and those derived here. For the TM sensors, the same 
parameters have been obtained and comparisons are made with the 
conventional bandwidth determinations. In addition, a spectral 
responsivity plot indicating the calculated parameters are 
presented for a typical TM channel . Since data from Band 6 
detector 4 was available at various temperatures, several temper- 
ature coefficients have been derived. 


RADIOMETRIC BANDWIDTH NORMAL I 2AT l ON 


Tha output signal (assumad a vol tags V) from a radiomatric sansor 
with a spactral rasponsivity ROJ whan axposad to a sourca 
(radianca from a surfaca, path radianca from an atmosphara, ate.) 
that producas a spactral irradianca E <x> at tha sanoor is 




ROJ dx 


< 1 ) 


In ganaral , wa wish to datarmina soma sourca charactaristlc , and 
cannot do this uni ass both tha sourca function and tha spac- 
tral rasponsivity ara Known ovar tha spactral intarval whara thoy 
ara both non-zaro. Howavar , if tha sp< ral rasponsivi ty R(X> 
can ba charactarizad as a ractangular function which has a valua 
Rn batwaan wavaMngth limits xl and X2 and 0 alsawhara, 
than wa may writa 


r 

V ■ Rn| E (X) 

J XI 


dx 


< 2 ) 


and tha sourca function is thus corractly avaluatad batwaan 

limits XI and X2. Radi oma trie bandwidth normalization is 

tha procass of assigning to a sansor with a Known spactral 


responsivity an equivalent squareband responsivity with well- 
defined wavelength limits and a constant responsivity Rn in the 
passband (Palmer, 1980). There have been at least 12 methods 
described throughout the literature, but only three will be 
considered here. 

The most widespread bandwidth normalization method, called peak 
normalization, arbitrarily sets the normalized responsivity Rn 
equal to the peak responsivity Rp of the sensor. The' bandwidth 
AK is then set such that the product Rn*4X is equal to 
the area under the spectral responsivity curve ROJ dx. 

This method tends to give results that are lower than true 
radiance levels and an ambiguity exists in that there is no well- 
defined means of locating the limiting wavelengths or defining 
the effective or center wavelength. 

Another popular method, called passband normalization, simply 
sets the wavelength limits at some arbitrarily-chosen respon- 
sivity level. A common level is the 50X point, and the resulting 
bandpass is called full-width at half-maximum < FWHM) . Since this 
method is primarily used only to define a bandwidth, a value for 
Rn is not normally assigned. If a center wavelength is given, it 
is usually the center of the passband, the average of the 
limiting wavelengths. This method was used by Markham and Barker 
(1983a, 1983b) . 
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The normalization technique used herein was -first described by 
Palmer and Tomasko ( 1980) and is called the moments method-. The 
analysis involves determination of the first and second moments 
of the spectral responsivi ty curve. The pertinent equations are: 

XC * | X*R(A)dA/J R (A ) dA 

<* 2 = [1 R(X)dxJ - A c * 

X2 * Xc ♦ jS-tf 
Xl * XC - 45‘0 

JX * X2 - Xl * 245-0 
Rn = (X2 - XI) -1 R( X) dX 

where 

xc = effective wavelength (centroid) 

02 = variance 

Note that this method gives fixed values for the wavelength 
limits, the effective (center) wavelength and the normalized 
responsivi ty with no ambiguities. The derivation, given in Palmer 
and Tomasko (1980), shows that exact results are achieved when 
the source function is quadratic, and superior results are 
obtained in comparison with other methods for continuous sources. 


okkjinal PAGE 19 

OF POOR QUALITY 




Figure 1 shows * typical moments normalization result, in this 
case -for TM -flight, band 2. Note that the e-f-fective wavelength is 
not located at the wavelength o-f peak response and that the limit 
wavelengths XI and X2 are not coincident with the 50% 
bandwidth points. 

I-f the spectral responsivity data available is relative rather 
than absolute, then the desired source -function cannot be deter- 
mined with Equation (2). Nevertheless, the bandpasses and effec- 
tive wavelengths calculated using the above method are still 
valid. This is indeed the case for the currently available MSS 
and TM data. When the absolute spectral responsivity at the 
wavelength of peak response is known, one need only to multiply 
it by the normalizing responsivity Rn to obtain the equivalent 
squareband response. Alternatively, the measured response to a 
known source can be used to determine Rn by integrating the 
source between wavelength limits XI and X2 and using 
Equation (2) . 
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RESULTS 


This section presents the results o-f the moments bandwidth 
normalization method as applied to the Landsat-4 and Landsat-D' 
MSS and TM scanners and compares them with previously derived 
values. For the MSS scanners, each detector was individually 
analyzed, and means and standard deviations were then determined 
-for each band. Table 1 shows the summary for the MSS Proto-flight 
instrument. The column headings are Rn for the equivalent square- 
band responsivi ty , Xl and X2 respectively for the short- 
and long- wavelength bandpass limits, dX for the bandpass and 
\c for the effective wavelength (centroid). The columns 
headed TM-83955 are the center wavelengths (calculated as the 
arithmetic mean of the 50X band limits) and the bandpasses (FWHM) 
as published by Markham and Barker (1983a). Table 2 gives the 
summary for the MSS Flight instrument. 

Comparison of the results given in Tables 1 and 2 show that the 
effective wavelength is the same for the two analysis methods 
within ±0.3X for the first three bands and about W. for band 
4. The effective passbands are somewhat different, with the 
moments analysis giving wider passbands than the 50X points 
(FWHM). The ratio of the passbands (FWHM/moments) is typically 
0.95 for the first three bands and 0.8 for band 4. 
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TABLE 1. MOMENTS SPECTRAL AfWLYSIS - MSS PROTOFLIGHT INSTRUMENT 


MOMENTS ANALYSIS i i TM-83935 


BAND 

: : Rn : 
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S JX 

: s xc s 
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TABLE 2. 

MOMENTS 

SPECTRAL ANALYSIS - MSS 

FLIGHT 

INSTRUMENT 




1 1 
1 1 

MOMENTS ANALYSIS 


S S TM-83955 S 


BAND 

! S Rn i 

XI S 

X2 

S Xc s 

4X 

S i Xc S 

JX 

1 

MEAN 

Si 91.18 S 

494.9 i 

610.7 

S 552.8 S 

116.0 

SS 551.8 S 

110.0 i 


C 

!! .55 ! 

1.35 S 

1.17 

i .37 S 

1.36 

SS .85 S 

.58 S 

2 

MEAN 

SS 91.14 S 

600.7 : 

698.8 

S 649.8 i 

98.4 

SS 649.8 S 

93.8 i 


<3 

SS .81 S 

.75 i 

.71 

i .44 i 

.77 

Si .42 S 
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3 

MEAN 

S ! 85.46 S 

698.5 S 
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S 756.7 ! 
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<J 

S : .63 S 

.61 : 

.25 

S .41 : 

.41 

s : 0 s 

0 i 

4 

MEAN 

S S 82.60 s 
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: 931.3 : 
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S! 922.4 ! 
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0 

S i 2.48 i 

1.9 S 

5.8 

! 3.7 ! 

3.7 

!! 5.7 i 

11.2 i 


J 
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‘V 


Table 3 shows the moments spectral analysis applied to the 
Proto* light TM instrument as f 1 own on Landsat— l. The thermal band 
(6) is treated separately and data are available for three 
temperatures for one detector (4). In this case, the 50X 
bandwidth points were determined by linear interpolation of the 
relative spectral responsivi ty tables in the vicinity of the 50X 
response points. The bandwidth is the FWHM and the effective 
wavelength xc is midway between the 50X points. Table 4 gives 
the summary results for the Flight TM instrument. 

Like the MSS data, there is little difference between the methods 
with regard to effective wavelength, but the moments method again 
gives wider passbands. Using data supplied by Barker ( 1982) , it 
was apparent that the detectors used for band 6 display strong 
temperature dependences , as shown in Figure 2. The relative 
spectral responsivi ty for Band 6, detector 4 is plotted at three 
different temperatures <90, 95 and 105K> . The temperature 

coefficients derived from this limited data set are shown in 
Table 5 with units in micrometers per Kelvin (except R) . 


- 8 - 


ORIGINAL PAGE « 
OF POOR QUALITY 


TABLE 3. MOMENTS SPECTRAL ANALYSIS - TM PR0T0FLI6HT INSTRUMENT 


6- 105K 


MOMENTS ANALYSIS 


450.3 

52 6.9 

621.3 

771.9 
1564 
2882 

10 . 2 ? 

10 . 2 ? 

10 . 2 ? 

18.32 

10.35 

10.22 


521.8 

615.6 

698.4 

906.3 

1791 

2351 


839.3 

1678 

2217 


71.5 

88.7 

77.1 

134.9 

227 

269 


1.680 

1.653 

1.676 

1.766 

1.566 


1507. BM ANALYSIS 


AX 


216.9 

250.2 


11.03 

11.04 
11.08 
11.18 
10.87 


NOTE: Band 6 data is for a temperature of 95K. Temperature data is 
for detector 4. 


TABLE 4. MOMENTS SPECTRAL ANALYSIS - TM FLIGHT INSTRUMENT 


MOMENTS ANALYSIS 


507. BN ANALYSIS 
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570.6 

660.7 
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82 
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! ! 11.43 1 
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:: 11.44 
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Table 5. TM Band 6 Temperature Coeff icients 


PARAMETER TEMPERATURE COEFF (/K) 


Rn ( responsivi ty> 

XI (short wavelength) 
X2 (long wavelength) 
A x (bandpass) 

Xc (centroid) 


- 1.12 
-.0086 F m 
-.021 
-.013 
-.014 


CONCLUSIONS 

An accurate radiometric bandwidth normalization method has been 
applied to the sensors on Landsat-4. It is recommended that the 
bandwidths and effective wavelengths presented herein be used in 
all cases where the source spectral radiance is unknown. It is 
further recommended that the moments method be employed to 
characterize future sensors. 
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Fig. 1. Moment* radiometric bandwidth normalization for a 
typical detector <TM PF Band 2) showing wavelength limits and Rn 
tor equivalent squareband response and effective wavelength. 


Fig. 2. TM Protoflight Band 6 , Detector *4 relative spectral 
responsivity curves at three temperatures. 





